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Abstract 

In colocated multiple-input multiple-output (MIMO) radar using compressive sensing (CS), a receive 
node compresses its received signal via a linear transformation, referred to as measurement matrix. The 
' samples are subsequently forwarded to a fusion center, where an -optimization problem is formulated 

and solved for target information. CS-based MIMO radar exploits the target sparsity in the angle-Doppler- 
range space and thus achieves the high localization performance of traditional MIMO radar but with 
0^ ■ many fewer measurements. The measurement matrix is vital for CS recovery performance. This paper 

' considers the design of measurement matrices that achieve an optimality criterion that depends on the 

• , coherence of the sensing matrix (CSM) and/or signal-to-interference ratio (SIR). The first approach 

(N 

. minimizes a performance penalty that is a linear combination of CSM and the inverse SIR. The second 

one imposes a structure on the measurement matrix and determines the parameters involved so that the 
SIR is enhanced. Depending on the transmit waveforms, the second approach can significantly improve 
\ SIR, while maintaining CSM comparable to that of the Gaussian random measurement matrix (GRMM). 

. Simulations indicate that the proposed measurement matrices can improve detection accuracy as compared 

to a GRMM. 
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I. Introduction 

Multiple-input multiple-output (MIMO) radar has received considerable recent attention |II]-|l3]. A 
MIMO radar consists of multiple transmit and receive antennas and is advantageous in two different 
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scenarios ||4]-||8l, namely, widely separated antennas and collocated antennas. In the first scenario ||4], 
the transmit antennas are located far apart from each other relative to their distance to the target. The 
MIMO radar system transmits independent probing signals from its antennas that follow independent 
paths, and thus each target return carries independent information about the target. Joint processing of 
the target returns results in diversity gain, which enables the MIMO radar to achieve high target resolution. 
Widely distributed MIMO radar systems are shown to offer considerable advantages for estimation of 
target parameters, such as location ||5] and velocity ||6]. In the collocated scenario QlE], the transmit 
and receive antennas are located close to each other relative to the target, so that all antennas view the 
same aspect of the target. In this scenario, the phase differences induced by transmit and receive antennas 
can be exploited to form a long virtual array with the number of elements equal to the product of the 
numbers of transmit and receive nodes. This enables the MIMO radar to achieve superior resolution in 
terms of direction of arrival (DOA) estimation and parameter identification Q. 

Compressive sensing (CS) theory Il9l- |[T2| states that a signal x that exhibits sparsity in some domain, 
can be recovered from a number of samples that is much smaller than that required by Nyquist theory. 
In particular, a signal of length that can be represented hy K (K « N) basis vectors in some space, 
can be recovered exactly with high probability from 0{K log N) measurements. Let *P denote the basis 
matrix that spans that space, and O denote an M x matrix with M ^ N, that is incoherent with 
T and is referred to as the measurement matrix. The recovery proceeds by finding the coefficients of 
the K basis vectors in the signal decomposition. This is formulated as an -optimization problem, i.e., 
min||s||i, s.t. to y = Ox = O'Ps. Throughput this paper, we will refer to recovery along these lines as 
the CS approach. The product O'P is usually referred to as the sensing matrix. According to the uniform 
uncertainty principle (UUP) ifTTI - llTSl . if every set of sensing matrix columns with cardinality less than 
the sparsity of the signal of interest is approximately orthogonal, then the sparse signal can be exactly 
recovered with high probability. In other words, CS recovery requires that O is incoherent with 'P. For 
an orthonormal basis matrix, use of a random measurement matrix leads to a sensing matrix that obeys 
the UUP with overwhelming probability |[TOl . The entries of such a measurement matrix can be taken 
from a Gaussian distribution or symmetric Bernoulli distribution. The rows of a Fourier matrix or an 
orthonormal matrix could also compose a measurement matrix. In this paper, we term as the conventional 
approach CS recovery using a Gaussian measurement matrix. 

The application of CS to radar and MIMO radar has been explored in |[14l-||T7l and |[T8l-ll211. 
respectively. In both |[T9l and ||20l . the authors considered a uniform linear array as a transmit and receive 
antenna configuration and proposed to use a submatrix of the identity matrix as the measurement matrix. 
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In ifTSl and ||2TI . a CS-based MIMO radar system implemented on a small scale network was proposed. 
The network consists of a number of transmit and receive nodes, each equipped with a single antenna, that 
are randomly distributed over a small area. Each transmit node transmits a different narrowband signal. 
If the number of targets is small, the signal that is reflected by targets and is picked up at a receive node 
is sparse in the angle-range-Doppler space. This fact can be exploited to achieve target detection and 
localization using only a small number of compressively obtained samples at each receive node, and/or 
by involving a small number of receive nodes lIlTI . The approach of IfTSl and ||2TI was applied to the 
case in which the targets are located within a small range bin and the sampling is synchronized with 
the first target return. To improve performance in the presence of strong interference the columns of the 
sensing matrix were designed to incorporate information on the transmit waveforms. 

In this paper, we consider a general scenario that does not confine the targets within a small range 
bin, nor does it require sampling synchronization. When the targets are separated by several range bins, 
different targets will introduce different delays in the received waveforms. In that case, the formulation of 
ifTSl and II2II no longer applies. This problem was considered in (221, where a step-frequency approach 
was proposed in order to improve range resolution. Here, our goal is optimal or suboptimal measurement 
matrix design that decreases the coherence of the sensing matrix (CSM) and/or enhances signal-to- 
interference ratio (SIR). The first design minimizes a performance penalty that is a linear combination 
of CSM and the inverse SIR. The measurement matrix is obtained by solving a convex optimization 
problem that involves high computational complexity. A suboptimal solution is also proposed that forces 
a specific structure to the measurement matrix. The second design targets only SIR improvement; it is 
constructed based on the transmit signal waveforms and accounts for all possible discretized delays of 
target returns within a given time window. It is shown that depending on the waveforms used, the latter 
measurement matrix can significantly improve SIR while it results in CSM comparable to that of the 
random Gaussian measurement matrix. 

The rest of the paper is organized as follows. In Section II we provide the signal model of a CS-based 
MIMO radar system with targets falling in different range bins. In Section III, we introduce the two 
proposed measurement matrices and provide the SIR analysis related to the second measurement matrix. 
Simulation results are given in Section IV for stationary targets. Finally, we make some concluding 
remarks in Section V. 

Notation: Lower case and capital letters in bold denote respectively vectors and matrices. The expecta- 
tion of a random variable is denoted by £"{•}. Superscripts i-)^ and Tr(-) denote respectively the Hermitian 
transpose and trace of a matrix. A{m, n) represents the (m, ?i)-th entry of the matrix A. Olxm denotes an 
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Lx M matrix with zero entries. 

II. Signal Model for CS-based MIMO Radar 

Let us consider a MIMO radar system consisting of M, transmit antennas (TXs) and A^^ receive antennas 
(RXs) that are randomly distributed over a small area (colocated). Each TX node transmits periodic 
narrowband pulses of duration Tp and pulse repetition interval (PRI) T. Let (H, ap/Crp aj) denote the 
location of the i-th transmit/receive node in polar coordinates. Let us also consider the presence of K 
slowly-moving point targets located in different range bins; the k-th target is at azimuth angle 6^ and 
moves with constant radial speed v^.. 

Let dk{t) denote the range of the k-th target at time t. Under the far-field assumption, i.e., dk{t) » r. , 
the distance between the ith transmit/receive node and the k-th target d'-Jd'-i^ can be approximated as 

4\t) « 4(0 - r]'/''(9k) = dm - rfl\eu) - vtt (1) 

where rj'/\dk) = r'/'^cos{dk - aj'). Let us consider the return from the ^-th target arriving at the /-th 
antenna during the m-th pulse, i.e., 

M, 

3^1(0 = Yj ^'^'^^ - (^'^-(^^ + expU27Tf{t - {d'.,{t) + d\,{t))lc)) (2) 

where c, / and yS/. denotes the speed of light, the carrier frequency, and the reflection coefficient of the 
^-th target, respectively; represents the transmit waveform of the i-th node. Under the narrowband 
assumption, and due to the slow target speed so that the Doppler shift is negligible, the baseband signal 
corresponding to ^ becomes 

M, 

yUt) « Yj ^^-^'^^ - 24(0)/c) exp(-727r/(4(0 + d\,{t))lc). (3) 
1=1 

Due to the closeness of the transmit and receive nodes, the distances between nodes and the target are 
approximately the same for all receivers. Thus, the time delay in the waveforms, induced by the k-th 
target can be approximately based on the range corresponding to the initial sampling time, i.e., dk{0), 
which is independent of the RX index. The /-th node compressively samples the return signal to obtain M 
samples per pulse (please refer to Fig. 1 of lIlTI for a schematic of the receiver). Let L denote the number 
of Tp/L-spaced samples of the transmitted waveforms within one pulse. The effect of the compressive 
receiver of Fig. 1 of ETI is equivalent to pre-multiplying by matrix O/ a T^/L-sampled version of the 
received pulse. The size of O/ is M x (L -i- L), where L is the maximum delay among the return signals 
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normalized by TpIL and is known in advance. Here M « L. The obtained samples are then placed in 
vector rim, which can be expressed in matrix form as ll22l 

K 

rim = A^'''"'"""^/D(A)Cr,Xv(et) + 0,n/„, (4) 

where 

1) Pimk - ~^''f^-^ + '^-^^+fk{m - l)T, where = ^ is the Doppler shift induced by the ^th target; 
X is an L X M, matrix that contains the transmit waveforms of antennas as its columns and 
diag{X«X} - [1,...,!]^; 

2) is the M X (L + L) measurement matrix for the /-th receive node; 

3) y{9k) - and D(/<;) = (liag{{e'^''f>^^'^"l^, . . _ ^e-i^^ML+l-DTpiL^y^ 

4) Tk = L cr /l J ^^'^ ^J"* ~ ^'^LxT^,iL,^Lx{L-Tt)^^ ■ Hcrc, we assume that the target returns completely 
fall within the sampling window of length (L + L)Tp/L, and that Tp/L is small enough so that the 
rounding error in the delay is small, i.e., x,(f - t/^) w Xi{t - ^jl)- 

5) n/„, is the interference of variance cr^ at the Z-th receiver during the m-th pulse, arising due to the 
jammer signals and thermal noise. 

Discretize the angle, speed and range space on a fine grid, i.e., respectively, . . . ,afjj, [hi,. . . ,bfjj 
and [ci, . . . ,cnJ. Let the grid points be arranged first angle-wise, then range-wise, and finally speed-wise 
to yield the grid points {an,b,„c„),n = I, ...,NaNbNc. Through this ordering, the grid point &„,^, c„ J 
is mapped to point (a„, bn, c„) with n = {tih - l)«a«c + ("c - ^)na + na- The discretized step is small enough 
so that each target falls on some angle-speed-range grid point. Then (01) can be rewritten as 

/ N 



rim = ^»/ 



Snei^'""-'\) Cl^jXv(«„) + ni„ 



(5) 



Pk, if the k-\h target is at (an,b„,c„) 
where N = NaNbNc, s„ = \ and 

0, otherwise 



-2c„f rj',{an)f 2bJ(m-l)T 

limn = + + ■ (d) 



In matrix form we have r;,„ = 0/mS -i- 0/11;^, where s = [si, ^a?]^ and 



Qim = O, [e'2'^'^""'D(2Z7i//c)C, Xy{cn), eJ^"'J""-D{2bNf/c)C, 2^ ,Xv(a^)] . (7) 



According to the CS formulation, 0/,„ is the sensing matrix and *P/,„ is the basis matrix. 



January 13, 2013 



DRAFT 



6 



If the number of targets is small as compared to A^, then the positions of the targets are sparse 
in the angle-speed-range space and s is a sparse vector. The locations of the non-zero elements of s 
provide information on target angle, speed and range. All the receive nodes forward their compressed 
measurements to a fusion center. We assume that the fusion center has the ability to separate the data of 
different nodes from each other. This can be done, for instance, if the nodes send their data over different 
carriers. The fusion center combines the compressively sampled signals due to Np pulses obtained at A^^ 
receive nodes to form the vector r. Using the predefined measurement matrices, the discretization of the 
angle-speed-range space, and also knowledge of the waveform matrix X, the fusion center obtains an 
estimate of s by applying the Dantzig selector |[23l . 



In this section, we discuss the design of the measurement matrix in order to improve the detection 
performance of CS-MIMO radar. For the sake of simplicity, we assume that all the nodes use the same 
measurement matrix, denoted by O, which does not vary with time. Since the targets are moving slowly, 
the Doppler shift within a pulse can be ignored. Generally, there are two factors that affect the performance 
of CS. The first one is the coherence of the sensing matrix. UUP requires low CSM to guarantee exact 
recovery of the sparse signal. The second factor is SIR. If the basis matrix obeys the UUP and the signal 
of interest s is sufficiently sparse, then the square estimation error of the Dantzig selector satisfies with 
very high probability Il23l 



where C is a constant. It can be easily seen from ([8]l that an increase in the interference power degrades 
the performance of the Dantzig selector. 

A. Measurement matrix design #1 

The goal of measurement matrix design is to reduce the coherence of the sensing matrix and at the 
same time increase SIR. The coherence of two columns of the sensing matrix, 0, corresponding to the 



III. Measurement matrix design 




(8) 
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k-th and k'-th grid point is given by 



)"tt'(0) 



PC 2c, Xv(a^) 



OC, 2,,, ,Xv(a,0 



I * I 



(9) 



where = C, 2^, XvCait). 

L77>7iJ 

Let the interference waveform at the l-th receive node during the m-th pulse be Gaussian distributed, 
i.e., «/„,(0 ~ CN{0,cr-). Let us also assume that the noise waveforms are independent across receive 
nodes and between pulses. Then the average power of the interference equals 



A',, N,. 



(10) 



m=l 1=1 



The average power of the echo reflected by the i-th target located on the fc,-th grid point of the angle -range 
space is approximately equal to 



(11) 



Therefore, the SIR equals approximately 



SIR 



(12) 



o-2Tr{0^0} 

The precise manner in which CSM and SIR affect the performance of the CS approach is unknown. 
Although theoretical bounds for the ^2 -norm of the estimation error have been proposed lE3l - ll25l . those 
bounds might not be relevant in applications in which the quantities of interests are the locations of the 
non-zero elements of the sparse signal, rather than the non-zero values themselves. This is the case in the 
problem at hand. In ll26l . an upper bound on the error probability of sparse support recovery, i.e., the total 
probability of missed detection and false alarm, under the optimal decision rule was derived. Although 
that upper bound is related to the detection of non-zero elements, it cannot be used for the design of the 
measurement matrix because it is rather loose, and further, it involves eigenvalues of submatrices of a 
given sensing matrix corresponding to all possible sparse patterns for the signal of interest. 

In this paper, we determine the measurement matrix by optimizing a linear combination of CSM and 
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the reciprocal of SIR. The CSM can be defined in various ways. Let us define CSM as the maximum 
coherence produced by a pair of cross columns in the sensing matrix. This criterion works well for a 
uniform sensing matrix but might not capture the behavior of the sensing matrix in cases in which the 
coherence of most column pairs is small ||27| . However, this coherence metric is widely used for the CS 
scenario due to its simplicity ||T7lll28l . 
The optimization problem becomes 



(13) 



where i is a positive weight, which reflects the tradeoff between the CSM and SIR. 

The problem of (fT2l) is not convex. In order to obtain a solution, let us first view (fT2T i as an optimization 
problem with respect to B = O^O. Furthermore, let us set the norm of the columns of the sensing matrix 
to 1, i.e., NrNpU^<&^<Duk =1,^=1, A^; this will significantly simplify the expression for iukk'{0) and 
Now, (fT2)) can be reformulated (using the approximation dEJ) as 



mm 



S.t. 



< t. 



t + iTr{B} 

m=l /=! 

k - \,...,N,k' = k+\,...,N 
N,NpU^Buk = l, k=l,...,N, 
B > 0, t>0. 



(14) 



which is a convex problem with respect to B. The first term in the objective refers to the maximum 
coherence of cross columns in the sensing matrix; the second term is proportional to the noise power 
which is a linear function of B. Once B is obtained, the solution of ([T2l ) can be obtained based on the 
eigendecomposition B = VLV^, as 



#1 



(15) 



where £ is a diagonal matrix that contains on its diagonal the nonzero eigenvalues of L, and V contains 
as its columns the corresponding eigenvectors. 

Another definition for CSM would be as the sum of the coherence of all pairs of columns in the 
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sensing matrix (SCSM). With this measure of CSM, the minimization problem becomes 



1 



mm 



SIR 



(16) 



Simulation results show that the solution of ([T6l ) can increase the number of column pairs in the sensing 
matrix that have low coherence as compared to the solution of (fT4)) . The latter solution tends to increase 
the coherence of some column pairs in the sensing matrix which had low coherence. This is intuitively 
expected because solving (fT4l ) requires more constraints than (fT6l ). In this paper, we use SCSM to design 
the measurement matrix. 

The above proposed methods for optimizing the measurement matrix reduce the coherence of cross 
columns in the sensing matrix without amplifying the interference. These methods improve the detection 
performance of the CS-based MIMO radar system, but incur an increased computational load as compared 
to CS-based MIMO radar that relies on the conventional measurement matrix. The number of complex 
variables entering the convex problem of ([16} is (^+^+^)(^+^) _ xhe computation complexity would be 
prohibitively high for large values of L + L. Also, for a large number of grid points A^, we have to deal 
with a large number of constraints. The optimal measurement matrix might be obtained and stored offline 
based on knowledge of grid points in the angle-range space. However, it would need to be updated once 
the basis matrix changes with the search area of interest. This would bring heavy burden to radar systems 
and thus might render the real-time application impossible. Therefore, ways to alleviate the computational 
load are worthy of investigation. A suboptimal scheme for this problem that involves lower complexity 
is discussed next. 

Let us impose a structure on the measurement matrix to be determined as follows: 



where W is an (L -i- L) X M unknown matrix to be determined and O is a M x MtiL + 1) Gaussian random 
matrix. Then the number of variables in W can be controlled by changing the value of M. We can obtain 
W by solving ([T6l ) with B = W^W and u^. = (DC, jc^ .X\{ak)- Furthermore, the structure in (ITTt enables 
a two-step processing for CS-based MIMO radar that simplifies the hardware of the receive nodes. In 
particular, a receive node linearly compresses the incident signal by using O. At the fusion center, all the 
signals forwarded by receive nodes are first multiplied by W and then jointly processed to extract target 
information. We can think of W as a type of post processing. In this way, the received nodes require no 
information about W, which reduces the communication overhead for the fusion center and nodes. 



0#i = WO 



(17) 
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In order to render the convex problem tractable, the norm of the columns in the sensing matrix is set 
to be a constant. This increases the number of constraints. If the number of variables is not sufficiently 
high, there might not be enough degrees of freedoms to decrease the coherence of the sensing matrix 
as compared to the original one. Since the number of constraints equals the number of grid points, the 
number of constraints can be decreased by reducing the search area. This can be done by considering 
grid points around some initial angle-range estimates, if such estimates are available. 

B. Measurement matrix design #2 

Although the suboptimal construction in ([TT] ) significantly reduces the number of variables, solving 
([T6l ) still requires high computational load. Further, the solution needs to be adapted to a particular basis 
matrix. In order to avoid these two shortcomings of 0#i , we next propose another measurement matrix 
that targets SIR improvement only. 

As in |[2n . we impose a special structure on the measurement matrix, i.e., 

0#2 = OW^ , (18) 

where O is an MxM (M < M) zero-mean Gaussian random matrix and W is an (L-i-L)xM deterministic 
matrix satisfying diaglW'^W} = [1,...,!]^. The above structure serves two purposes. First, the matrix 
W can be selected to improve the detection performance of the CS approach at the receiver. Second, 0#2 
is always Gaussian regardless of W. As will be shown next, with the appropriate W, 0#2 may not result 
is higher CSM as compared to the conventional measurement matrix. Next, we discuss the selection of 
W. 

The average power of the echo reflected by the ^-th target with respect to node locations, conditioned 
on the transmit waveforms, is approximately equal to 

P', * tS*pA^,A^^,£{Tr{OW^C,,Xv(6^,)(OW^C,,Xv(0,))^}} 

w t6*pA^,A^/,Tr{OW^C^,,X£{VA.}X^C^WO^} (19) 

where V^. = v(%)v^(0^) and its {i,j)-t\\ entry can be expressed as Vk{i,j) - 

As already noted, the Doppler shift within a pulse is ignored in ([T9l ). Since the nodes are uniformly 
dispersed on a disk of radius r, the distribution of h = -y sin(a. - i/'o) is given by |[29l 

fh{h) = - Vl - h^, -l<h<l (20) 



January 13, 2013 



DRAFT 



11 



so that 



E \e'"^'\ = 2 



a 



(21) 



where is the first-order Bessel function of the first kind. Thus, based on (|2T] ) we can obtain |[2TI 



^(4sin(^)) 



= j and k = k' 
= j and k k' 



(22) 



where ^(jf) = 2 ' "^.^ . As observed in ETI . the terms multiplied by $"^(2) are small enough and can be 
neglected. Thus, the average power P\ in ([T9l ) can be further approximated by 



-TrlW^C^XX'^C^W}. 



^.pMiV,iV,^__,„ „^ ^^^^^^ 



(23) 



M 

Inserting <D#2 into (fTOl) . the average power of the interference can be approximated as 

M 

P„ - cr^NpN,.TmW"W^"] = cr^NpN,- 3>g/w,-,0;. * a^NpN^M (24) 

?=i '.,/■ 

where O,, and w,, are the (/, 7)-th entries of and W^W, respectively. The approximation in (|24] | uses 
the constraint diag{W"W} = [1, . . . , 1]^ and the fact that Iji^j^qiWij*^*^j ~ for sufficiently large 
M due to Og, ~ AA(0, 1/M). 

Based on (l23l) and ((24l) . the SIR is given approximately by 

|2 



5/7?, - « ^Tr{W^Q,,W) 



(25) 



where Q^^ = Cr^XX^C^ is an (L + L) x (L + L) matrix of rank M,. The maximization of SIRk over W 



can thus be approximated by the problem 



W* = max ^^Tr{W"Q^,W} 

W,M cr^M 

s.t. diag{W^W} = [l,...,l]^^j. 



(26) 



It can be easily seen that W* contains as its columns the eigenvectors corresponding to the largest 
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eigenvalue of Q^^j. Since the largest eigenvalue of Q^^^ is not greater than TrlQ^J = M,, the maximum 
SIRk is bounded by 

Boun,l:^,SIR,,^. (27) 

The upper bound is achieved when the rank of X equals 1, i.e., all the transmit nodes send out the same 
waveforms. When orthogonal waveforms are utilized, i.e., X^X = Im,, the SIR^ reaches the lower bound. 

It can be shown that, when the transmit waveform are orthogonal, i.e., X^X = 1^;, Qrj. has M, nonzero 
eigenvalues which are all equal to 1. Therefore, for a fixed M, M < M„ the optimal W contains the M 
eigenvectors of Qk corresponding to eigenvalue 1 and achieves maximum SIRk equal to SIRk = 
Since the maximum 5 IRk is independent of M, any matrix containing M,M < M,, eigenvectors of Q 
corresponding to eigenvalue 1 would give rise to the maximum 5/7?^.. However, M - M, results in smaller 
CSM than any M less than M, due to the fact that the rank of W is M. Therefore, the optimal W is 

W" - C^,X. (28) 

For the case of completely coherent transmit waveforms in which the upper bound in (l27l ) is achieved, 
the resulting W** is rank deficient. 

Unfortunately, W** is not achievable since it depends on the time delay induced by a target, Tk, which is 
unknown. To address this issue, we replace SIRk in the objective function in (l26l ) with the average SIRk, 



where the average is taken over all possible delays, and is denoted here by SIRk- Assuming that the time 
delay induced by the k-t\\ target follows a discrete uniform distribution, i.e., p{Tk - k) - J-r, k = 0, . . . ,L, 



we can wnte 



SIRk = y ^^Tr{W^Q,W) - ^^'L , ^ Tr{W^CW} (29) 

cr^M ^ L + 1 cr^M L + 1 

r=0 



where 



C = Yj^r = [CoX, . . . , QX][CoX, . . . , CiXf. (30) 



T=0 



Therefore, the optimization problem that maximizes SIRk can be rewritten as 



W* - m&xSIRk 

W,M 



diag{W^W} = [l,...,l]^^^. (31) 
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The solution W* of the above problem contains as its columns the eigenvectors corresponding to the 
largest eigenvalue of C. Unlike (l26l ). we cannot find a close-form solution to (OTI ) that has sufficiently 
high rank. Further, the problem of (|3Ti is non-convex. Inspired by the form of (|28T l. we propose a feasible 
W by taking all possible delays into account as follows: 

W = [CoX,...,CiX]. (32) 

Since C,X contains eigenvectors corresponding to the largest eigenvalues of Q,-, utilizing (l32l ) results in 
the average S IR^ bounded by 

Bound!: — tt- i-A<5/«a< — + A (33) 

0-2 L 1 0-2 L 1 

where A denotes ^.fif ,., Tr{y.,,^,X^C^Q^.C^X}. One can see that Bound 2 would be reduced to Bound 
1 when L = 0. 

Next, we examine the resulting SIR based on three types of waveforms, namely a rectangular pulse, 
independently generated quadrature phase shift keying (QPSK) waveforms and Hadamard codes. In 
particular, we show that using 0#2 can suppress interference uncorrelated with the transmit waveforms, 
and maintains coherence as low as that corresponding to the Gaussian random measurement matrix. 

I) SIR under the conventional measurement matrix: Let us consider a conventional measurement 
matrix O^, which is an M x (L -i- L) Gaussian random matrix of unit column norm with TrjO^Of } = M. 
The average power of the interference is Pn = cr^M (see (fTOll). 

Let S, be a square matrix, formed by shifting the main diagonal of 1^ up by / positions. It can be easily 
seen that = S_/. The average power of the target returns from K targets at a receive node, conditioned 
on the transmit waveforms, equals 

P, = E{r"r\X} - 2 + Z ^^"^^ 

k ki^k' 
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where 



P\ = ^,p£{Tr{O,C,,Xv(%)(OC,,Xv(0,))^}} 



|2 



L + L 



(35) 



and 



Pf'^ «^*^,,^2^4sin(^^)). __Tr{X^S,,_,,X}. (36) 

" s, ' 

Ttt' 

2) SIR for the measurement matrix 0#2.' The proposed measurement matrix 0#2 - ^>W'^ results in the 
same average interference power as the matrix O^,. The average power of the desired signal conditioned 
on the transmit waveforms, P,, however, will improve. Like (|34] |. can be partitioned into the sum of 
the autocorrelation, P\, and cross correlation, Pj'*^ , of the returns from K targets. It holds that 

P\ = ^,p£{Tr{0#2C,,Xv(ei)(6#2C,,Xv(e,))«}} 
w ^;tpTr{OW^C^,XX"C^WO^} 

^ TrlW'^Cr.XX^C^^W} 



(L + l)M, 
and 

7kk'M ^ 



y Tr{X^S„_^,XX^Sf ^ X} (37) 



pkX ^ Jkk'M y xr{x^S,,_,XX^S,_,,,X}. (38) 
(L + 1)M, Aj ' 1 ^ I 

For orthogonal, or randomly generated waveforms across the transmit nodes, P^ always dominates the 
average power of the desired signal. In order to increase P^, the quantity TrlX^S^-r^XX^S^^^X} in (l37l) 
needs to be as large as possible. X^S„,X can be expressed as 

i X" X„,+i.L, m > 
X«S,„X= (39) 

[ Xf_,^,.^Xi:z,+„„ Otherwise 

where X,; , denotes the matrix that contains the rows of X indexed from / to / 

Eq. (l39l ) implies that the non-circular autocorrelation of the waveform sequence of a transmit node, 
i.e., Riir) = jj^ Xi{t)x*.{t - T),i = 1, . . . , M,, should be insensitive to the shift. This essentially requires a 
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narrowband signal. Based on this principle, the best candidate is a rectangular pulse and the maximum 
Ps equals 

^ ^ MM, y y |. ,2 ( L-\q-Tt\ \^ ^ M_M, ^^{L- |t, - q\){L -\q- t,,\) 



( K \ 

<MM, Yj^^kf + T^ykk' 



(40) 

a=l ktk' 

The equality in ( |40l ) holds only if the targets induce identical delays and the sampling window is of a 
length that exactly covers the duration of target returns. Obviously, the transmit nodes cannot use identical 
waveforms. This is because the transmit waveforms are required to be orthogonal, or randomly generated 
in order to maintain low CSM. 

Similarly, the minimum average power of the desired signal is achieved when randomly generated 
QPSK waveforms are used, because such waveforms cover the widest bandwidth for the fixed pulsed 
duration Tp and the length of waveforms L. The corresponding value of Ps is approximately 



(L+l)fY 

^ ■' k=\ \q=0,qi^Tt 



(41) 



For orthogonal Hadamard waveforms that are of smaller bandwidth than the randomly generated QPSK 
waveforms, the average power of the desired signal equals approximately 



■7= 

Recall that W* corresponding to the true delay gives rise to the maximum received signal power. Adding 
the terms Cf^X, f;^ ?t r^, to W (see (l32l )) would lower P,. When a coarse delay estimate is available, we 
need to consider only delays around that estimate and thus the length of the sampling window can be 
shortened. This effectively reduces the number of possible delays that are considered for the construction 
of W. Therefore, can be improved for the waveforms considered above if a coarse delay estimate is 
available. 

3) The SIR gain: Let SIRp and SIRc denote the SIR corresponding to measurement matrices 0#2 
and Oc, respectively. When transmitting Hadamard codes, the SIR gain induced by using the proposed 



January 13, 2013 



DRAFT 



16 



measurement matrix can be expressed as 



M x^K x^L \o a i L-\q-Tt\ \^ 
2^k=l L+L 

^ L + L Zf^i lAf q 



(43) 



where 



Q- = J] (L - 1^ - T,|)2 = (L + 1 - 2L)(T, - 1/2)2 + ^ ^2 + ^ . (44) 

9=0 g=L-L 

For a fixed Z and L, with < L < 2L - 1, Q can be bounded as 

t^'^Ck< j^,^ + ^2^-l_^. (45) 

q=L—L q=L-L 

Therefore, lower and upper bounds on the approximate SIR gain using Hadamard codes are given by 

~ ' , < THa < ^ ^ = . (46) 

MXL + 1)L2 M,(L + 1)L2 

Similarly, the SIR gain using randomly generated QPSK waveforms is bounded by 

{L + z:)(i:^^^_, ? + f - L) (L + L)(i:^^,_, ^ + f + f - l) 

^ ■ < rnpsK < — z ■ ■ (47) 

(L+1)L2 (L+1)L2 

Z"- - q- 

As long as L < L and M, < L, vqpsk is always greater than 1. When i ^ — > 1, the lower bound on rna 
is higher than that on tqpsk- For a sufficiently long L and moderate M„ r^a would be superior to rgpsj^. 
Based on (l40l ) and (l42l ). one can infer that the SIR gain using the rectangular pulse is approximately Mr 
times greater than that using Hadamard codes. 

4 ) The CSM based on the suboptimal measurement matrix #2 : In this section, we examine the effect 
of the proposed W in (l32l ) on the CSM. For simplicity, the targets are considered to be stationary and 
the possible delays for constructing W are based on the range grid points used to form the basis matrix. 
Then the sensing matrix based on 0#2, or the Gaussian random matrix can be respectively represented 
as 

= O^^'P = W^WV (48) 
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and 

= ^MXCL+D^P = ^MX(L.I)WV (49) 

where V = kron{Ii^i, [\{ai), . . . ,\(aNj]) and 0,x/ is an / x j Gaussian random matrix whose entries are 
of zero mean and variance 1 /j. For sufficiently large j, the column coherence of can be approximated 

as 

m'i®) ^ I (50) 

V2Ki:«, "^(m, iW(m, i))vuk{i) Z/d]™ 0(m, /)0*(m, i))vu'k'ii) 
where v^'CO denotes the /-th diagonal element of the matrix W^C, ic^ |Xvn;(a^.)| W^C, 2.,., ,Xvm(a^,0| . 

'-'■rp/L-l \ ^ 

Without loss of generality, we let the columns of O be of unit norm. Then (|50] l can be further written as 

m'(0) = , (51) 

One can easily see from (ISTI ) that the coherence of is approximately equal to that of matrix W^WV. 
The same conclusion applies to as well, i.e., the coherence of is approximately equal to that of 
matrix WV. Since W^W is more ill-conditioned than W, the conditional number of W^WV is greater 
than that of W^W. Therefore, using 0#2 increases the maximum CSM as compared to the Gaussian 
random measurement matrix with high probability. However, for a well conditioned W, the increase of 
the maximum CSM caused by 0#2 is negligible. 

C. 0#i v.s. 0#2 

We have proposed two measurement matrices based on different performance metrics. The advantages 
and disadvantages of 0#i and 0#2 are summarized as follows. 

• Complexity 

Solving 0#i involves a complex optimization problem and depends on a particular basis matrix, while 
0#2 requires knowledge only of all the possible discretized time delays. Therefore, the construction 
of 0#i involves higher computational complexity than does 0#2. 

• Performance 

0#i aims at decreasing the coherence of the sensing matrix and enhancing SIR simultaneously. 
The tradeoff between CSM and SIR results in 0#i yielding lower SIR than 0#2. Therefore, 0#i is 
expected to perform better than 0#2 in the case of low interference, while it should perform worse 
in the presence of strong interference. 
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IV. Simulation Results 

In this section, we demonstrate the performance of CS-based MIMO radar when using the proposed 
measurement matrices 0#i and 0#2, respectively. We consider a MIMO radar system with transmit and 
receive nodes uniformly located on a disk of radius 10m. The carrier frequency is / = 5GHz. The 
received signal is corrupted by zero-mean Gaussian noise. The signal-to-noise ratio (SNR) is defined as 
the inverse of the power of thermal noise at a receive node. A jammer is located at angle 7° and transmits 
an unknown Gaussian random waveform. The targets are assumed to fall on the grid points. 

A. The proposed measurement matrix 0#2 

1) SIR improvement: M = 30 compressed measurements are forwarded to the fusion center by 
each receive node. The maximum possible delay is L = 100. Figure \T\ compares the numerical and 
theoretical SIR produced using the rectangular-pulse, Hadamard waveforms and randomly generated 
QPSK waveforms for the case of M, = 30 transmit nodes and Nr = I receive node. The SIR performance, 
shown in Fig. [1] is the average of 1000 independent and random runs. The theoretical SIR of these three 
sequences is calculated based on ( |40l ). (I42l ) and (|4T]) . respectively. The power of thermal noise is fixed 
to 1 and the power of the jammer varies from -20dB to 60dB. Applying the proposed measurement 
0#2 at the receivers, the rectangular pulse and Hadamard waveforms produce a significant SIR gain 
over the Gaussian random measurement matrix (GRMM), while the random QPSK sequence achieves 
almost no gain. Furthermore, the numerical SIR performance follows the theoretical SIR for all three 
sequences. Figure |2] demonstrates the SIR performance obtained by averaging over 500 independent runs 
corresponding to independent interference waveforms, for diff'erent values of the maximum time delay 
L. We consider a case in which only one target exists and the jammer power is 225. One can see that 
a decrease in L can significantly improve SIR yielded by QPSK waveforms when L is less than 10. In 
contrast, Hadamard waveforms and rectangular pulse produce almost the same SIR for different values 
of L. This indicates that for QPSK waveforms the prior information of possible delays enables SIR 
improvement, while for the other two types of waveforms considered prior information did not make a 
difference. 

2) The CSM: Figure |3] shows the histograms of the condition number and the maximum CSM using 0#2 
for Hadamard waveforms and the GRMM produced in 100 random and independent runs. We consider the 
case of M = 30, M, = Nr = 10 and the grid step of the discretized angle-range space is [0.5°, 15m]. One 
can see that the numerical results fit the observations in Section IIlI-B4[ i.e., 0#2 increases the maximum 
CSM as compared to the GRMM with high probability. In Fig. |4] we use histograms to compare the CSM 
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corresponding to adjacent columns over 100 independent and random runs. Although 0#2 incorporates 
information about the waveforms, the distribution of the column correlation does not change significantly 
as compared to that of the conventional matrix. Among the three types of waveforms, the rectangular 
pulse gives rise to the worst CSM distribution, indicating that the performance of the proposed CS 
approach would be significantly degraded if rectangular pulses are transmitted. This is because the high 
autocorrelation of the rectangular pulse results in high CSM independently of the measurement matrix 
used. 

B. The proposed measurement matrix 0#i 

We consider a scenario in which M, = N,- = 4- and three stationary targets exist. The azimuth angle and 
range of three targets are randomly generated in 100 runs within [0°, 1°] and [1000m, 1090m], respectively. 
The data of only one pulse is used and thus only the angle -range estimates can be obtained. The spacing 
of adjacent angle-range grid points is [0.2°, 15m]. 0#i is obtained from (IT4l ) based on the special structure 
of (fTTl) . O in ([TT] ) is replaced with 0#2. We consider diff'erent values of the tradeoff" coefficient A in ([T4l ). 
Transmit nodes send Hadamard waveforms of length L = 128. Only M = 20 measurements per pulse 
are collected and forwarded to the fusion center by each node for CS-based MIMO radar while 100 
measurements are used by the MIMO radar based on the matched filter method (MFM) Il30l . 

Figure |5] shows the distribution of CSM for the GRMM, <D#i and ^>#2 in 100 random and independent 
runs. One can see that the GRMM and C»#2 lead to similar coherence distributions. 0#i slightly reduces the 
maximum CSM and significantly increases the number of column pairs with low coherence as compared 
to the other two measurement matrices. 0#i obtained from (IT6] l using A = 0.6 and ^ = 1.5 produce a 
similar coherence distribution. Figure |6] shows the SIR performance of CS-based MIMO radar using the 
GRMM, 0#i and 0#2, for diff'erent values of noise power in the absence of a jammer. One can see from 
Fig. [6] that 0#2 outperforms the other two measurement matrices in terms of S IR. 0#i obtained from 
([T6l ) using A = 0.6 yields slightly better SIR than GRMM. As expected, increasing A from 0.6 to 1.5 
moderately improves SIR. 

Figures 17] and [8] compare the ROC performance of CS-based MIMO radar using the three aforemen- 
tioned measurement matrices and MIMO radar based on the MFM, for diff'erent combinations of SNR 
and jammer-signal power. The probability of detection (PD) here denotes the percentage of cases in which 
all the targets are detected. The percentage of cases in which false targets are detected is denoted by the 
probability of false alarm (PFA). It is demonstrated in Figs. |7] and [8] that 0#i and 0#2 with Hadamard 
waveforms can improve detection accuracy as compared to the GRMM in the case of mild and strong 
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interference, respectively. Since an increase in the tradeoff coefficient A can enhance SIR, 0#i obtained 
from ([T6l ) using 1 = 1.5 performs better in the case of strong interference than using A = 0.6. Note that 
the three measurement matrices give rise to similar performance for SNR = lOdB and yS = 0. This is 
because the interference is sufficiently small so that all the measurement matrices perform well. Again, 
one can see that the MFM is inferior to the CS approach although it uses far more measurements than 
the CS approach. 

It has been seen from Figs. |7] and [8] that the tradeoff coefficient A affects the performance of CS- 
based MIMO radar using 0#i. In order to further investigate the effect of A, the curves of probability 
of detection accuracy are shown in Fig. |9]for 1 = 0.6, 1, 1.5, 2 for different thresholds of hard detection. 
The probability of detection accuracy here denotes the percentage of cases in which no real targets are 
missing and no false targets exist. By taking all four combinations of SNR and jammer-signal power into 
account, 1 = 1.5 results in the best performance. For a particular case, the optimal tradeoff coefficient 
depends on multiple factors, i.e., the basis matrix and the interference. The manner in which SIR and the 
CSM affect the support recovery of a sparse signal still remains unknown. Therefore, it is impossible to 
theoretically determine the optimal tradeoff coefficient. 

V. Conclusions 

We have proposed two measurement matrices in order to improve target detection performance of 
CS-based MIMO radar for the case in which the targets may be located across several range bins. 
The first one 0#i aims at enhancing SIR and reducing the CSM at the same time. It is obtained by 
solving a convex optimization problem. This measurement matrix requires a heavy computational load 
as compared to the conventional measurement matrix, and also needs to adapt to a particular basis 
matrix. The computational burden of solving 0#i can be alleviated through reducing the number of 
variables involved in the optimization problem. The second proposed measurement matrix 0#2 targets 
improving SIR only. It is constructed based on the transmit waveforms and also accounts for all possible 
discretized delays of target returns within the given time window. ^>#2 is dependent on the range grid only 
and requires much lower complexity than does It is shown that 0#2 based on reduced bandwidth 
transmit waveforms can improve SIR, but on the other hand, using waveforms that are too narrowband 
increases the CSM, thus invalidating conditions for the application of the CS approach. Therefore, the 
waveforms must be chosen carefully to guarantee the desired performance using the second measurement 
matrix. Numerical results show that ^>#i and 0#2 with the proper waveforms (e.g., Hadamard codes) can 
improve detection accuracy as compared to the Gaussian random measurement matrix in the case of 
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small and strong interference, respectively. 
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Fig. 1. SIR corresponding to GRMM and <D#2 for different transmit waveforms (M = M, = 30 and A', = 1). 
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Fig. 2. SIR corresponding to <I>#2 for different values of L (M = M, = 30 and = 1). 
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Fig. 3. Conditional number and the maximum coherence of the sensing matrix based on 0#2 (M = 30 and = M, = 10). 
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Fig. 4. Coherence of adjacent columns of the sensing matrix based on <D#2 for different transmit sequences (M = 30 and 
M, = N, = 10). 



January 13, 2013 



DRAFT 



25 



GRMM 

2000 1 , , , , , — 



1000 



lb 



1 1^ 



ILii 



1000 
0, 



ii< ijLii 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

$„, tradeoff coefficienl=0.6 

#1 

2000 
« 1000 

E °0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
2 $„, tradeoff coefficient=1 .5 

O) " 

B 2000 1 1 1 1 1 1 1 1 1 1 1 



1000 

% 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
2000 1 1 1 1 1 1 1 1 1 1 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

CSM 



Fig. 5. Coherence distribution of cross columns of tlie sensing matrix using <D#i, <D#2 and the GRMM (M, = A', = 4). 
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Fig. 6. SIR for CS-based MIMO radar using 0#i, <I>#2 and GRMM for different values of noise power (M, = A', = 4). 
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Fig. 7. ROC curves for CS-based MIMO radar using <D#i, 0#2 and the GRMM and for MIMO radar using the MFM 
(M, = yV,. = 4 and A = 0.6). 
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Fig. 8. ROC curves for CS-based MIMO radar using <D#i, <D#2 and the GRMM and for MIMO radar using the MFM 
(M, =N, = 4andA= 1.5). 
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Fig. 9. The probability of detection accuracy for CS-based MIMO radar using 0#i, 0#2 and the GRMM for different values 
of A (M, =Nr = 4). 
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